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Abstract 

In the present paper we construct stock price processes with the same marginal log- 
normal law as that of a geometric Brownian motion and also with the same transition 
density (and returns' distributions) between any two instants in a given discrete-time grid. 

We then illustrate how option prices based on such processes differ from Black and 
Scholes', in that option prices can be either arbitrarily close to the option intrinsic value 
or arbitrarily close to the underlying stock price. 

We also explain that this is due to the particular way one models the stock-price 
process in between the grid time instants which are relevant for trading. 
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The theoretical result concerning scalar stochastic differential equations with pre- 
scribed diffusion coefficient whose densities evolve in a prescribed exponential family, 
on which part of the paper is based, is presented in detail. 

Keywords 
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1 Introduction 

Since the seminal papers by Black and Scholes (1973) and Merton (1973), many researchers have 
tackled the problem of option pricing and hedging by resorting to continuous-time mathematics 
to fully exploit the richness of its theoretical results. In particular, the underlying stock price 
dynamics has often been modelled through a geometric Brownian motion. This assumption is 
still quite popular, especially among practitioners, even though the empirical distribution of 
stock returns is often found more leptokurtic than the normal one. Indeed, the tradeoff between 
analytical tractability and empirical findings still makes the geometric Brownian motion a quite 
acceptable choice. 

However, as we shall prove in this paper, the geometric Brownian motion is not the only 
continuous-time process that possesses a lognormal marginal distribution and normal log- 
returns. In fact, we shall construct a family of processes with these characteristics and then 
we shall also analyse the main implications of the existence of such a family as far as option 
pricing is concerned. 

In order to achieve our goal, we consider scalar nonlinear diffusion processes whose densities 
evolve in given finite-dimensional exponential families, as from Chapter 7 of Brigo (1996). We 
begin by a stochastic differential equation (SDE) a la Black and Scholes for the stock price, 
whose solution is the traditional geometric Brownian motion. The probability density of such a 
stock price model evolves in the exponential family of log-normal densities. It is then possible 
to define a second SDE with a different (and arbitrarily given) diffusion coefficient such that 
its solution has the same marginal lognormal density as that of the Black and Scholes (1973) 
geometric Brownian motion. 

Furthermore, one can choose the drift of this second SDE in such a way that its solution has 
transition densities, between two any instants of a prescribed finite subset of the time interval, 
that match the transition densities of the Black and Scholes process. We thus derive a family 
of stock price processes that behave almost cquivalently to a geometric Brownian motion in 
that they are probabilistically identical along the finite subset of the time interval. 

We then consider the problem of option pricing in a continuous-time framework. The 
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interesting result we obtain is that notwithstanding the above mentioned affinity with the 
Black and Scholes model, the stock price processes of our family lead to completely arbitrary 
option prices, ranging the whole interval of prices inbetween the no-arbitrage lower and upper 
bounds. Such a seeming paradox can be explained in terms of the differences between the 
dynamics of these stock price processes at an infinitesimal level. 

Our result is similar in spirit to that of Rogers and Satchell (1996). However, our approach 
is much more constructive in that we are able to provide, under the original measure, easy 
and explicit dynamics for the stock price that are theoretically consistent with whatever option 
price is observed in reality. 

The practical implications of our result can be analysed by considering the case of a prac- 
titioner who has to price an option. Although the physical time he works in is discrete, he 
usually resorts to continuous-time mathematics modelling the stock price underlying the op- 
tion through a Black and Scholes (1973) process. However, as the previous theoretical result 
implies, the geometric Brownian motion is only one of infinitely many processes that possess 
the discrete-time properties required by the practitioner. Therefore, as far as the stock price 
dynamics is concerned, the practitioner must view all these processes as equivalent to each 
other. Yet, when option pricing is taken into account, he must also be aware that they all 
imply different, and arbitrary, option prices. Since this paradoxical situation reveals a clear 
limit of the approximation of discrete time with continuous time, our practitioner can react as 
follows. He simply avoids the theory of market completeness and just chooses a process which 
is consistent with the option price he believes in. Such an option price can be either exoge- 
nously given by the market or endogenously produced by a model for pricing and hedging in 
incomplete markets. This is a further contribution of the paper. We do not derive option prices 
and hedging strategies assuming a specific evolution of the underlying stock price. Instead, we 
define a family of processes with equivalent characteristics among which one can select a single 
process according to his pricing and hedging purposes. 

The paper is structured as follows. In Section 2, the theoretical result on the construction 
of some suitable diffusion processes is briefiy reviewed. Section 3 deals with the application of 
this result to the case of the Black and Scholes (1973) model. Section 4 analyses the problem 
of option pricing in continuous time and considers a natural particular case. Section 5 analyses 
some practical implications of using any of our alternative processes in describing the asset 
price dynamics. Section 6 concludes the paper. 

A reduced version of this paper can be found in Brigo and Mercurio (2000). 
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2 Stochastic differential equations and exponential fam- 



In this section we consider the technical result the paper is based on. This result solves the 
following problem for scalar diffusion processes: Given a diffusion coefficient and a curve in an 
exponential family of densities, find a drift such that the solution of the resulting stochastic 
differential equation (SDE) has a density evolving in the prescribed exponential family according 
to the given curve. This problem has a straightforward solution. In this section, we shall present 
a short summary of the steps leading to the problem formulation and to its solution. This 
summary is based on Chapter 7 of Brigo (1996). A related result appeared in Brigo (1997), 
and the general result with its applications both to mathematical finance and to stochastic 
nonlinear filtering is also reported in Brigo (2000). 
Let us consider the SDE 



where {Wt,t > 0} is a standard Brownian motion independent of the initial condition Xq. 
In order to make sure that we are dealing with an equation whose solution exists unique, we 
formulate the following assumption. 

(Al) The initial condition Xq is a continuous random variable with density po{x) > for all 
X G R w.r.t. the Lebesgue measure on R whose moments of any order are finite. Moreover, 
the stochastic differential equation ([1]) characterized by the coefficients /, a, and by the 
initial condition Xq admits a unique strong solution. 

Explicit conditions ensuring (Al) are, for example, local Lipschitz continuity and linear 
growth, or the Yamada-Watanabe condition (see e.g. Rogers and Williams (1987), Section 



Once existence and uniqueness of the solution of a SDE have been established, we can 
analyse the distribution of its solution at all time instants. In describing the evolution of 
the distribution of a diffusion process, the Fokker-Planck partial differential equation is a 
fundamental tool. We therefore introduce the following assumption. 

(A2) The unique solution Xt of ([H) admits a density pt that is absolutely continuous with 
respect to the Lebesgue measure, i.e.. 



ilies 



dX, = ft{Xt)dt + at{Xt)dWt 



(1) 



V-40). 




A 



Pt{x)dx, for all Borel sets A, 



and that satisfies the Fokker-Planck equation: 



dpt 
dt 



d_ 

dx 



iftPt) + 




(atpt), at{-) = a^i-) . 



2 Q^2 



D. Brigo, F. Mercurio. Banca IMI - PDG internal report 



5 



Examples of assumptions on the coefficients /, a and on their partial derivatives under which 
(A2) holds are given in the literature. See for example Stroock and Varadhan (1979) or Fried- 
man (1975). 

In order to appropriately introduce the problem we mentioned at the beginning of the 
section, we now present a definition of exponential family. 

Definition 2.1 Let {ci, ■ ■ ■ , Cm} be scalar functions defined on R, such that {1, ci, ■ ■ ■ , Cm} are 
linearly independent, have at most polynomial growth, are twice continuously differentiable and 
the convex set 

00 := 1^ = . . . , 6""}' G R"" : ^{6) = log J exp[0'c(x)] c/x < ooj , 
has non-empty interior, where c{x) = {ci(x), ■ ■ ■ , Cmix)}' and " denotes transposition. Then 

EM{c) = {p{-, 9),ee 0}, p(x, 6) ■= exp[^'c(x) - V'(^)] , 
where B C Bq «s open, is called an exponential family of probability densities. 

Our problem consists in finding a SDE whose solution Xt has a density pt that follows a 
prescribed curve in a given exponential family. More precisely, we require the curve t ^ pt, in 
the space of all densities, to coincide with a given curve t ^ p{-, 6t) in a given i?M(c)ll] 

This is formalized in the following. 

Problem 2.2 Let be given an exponential family EM{c), an initial density pq contained in 
EM{c), and a diffusion coefficient at{-) := cr^(-); ^ > 0. LetlA{pQ,a) denote the set of all drifts 
f such that Pq, f and a and its related SDE (Q]) satisfy assumptions (Al) and (A2). Assume 
U{po,cr) to be non-empty. 

Then, given the curve t i— > p{-,Ot) in EM{c) (where t ^ 6t is a C^-curve in the parameter 
space B ), find a drift in U{pq, a) whose related SDE has a solution with density pt = p{-, 6t). 

The solution of this problem is given by the following. 

Theorem 2.3 (Solution of Problem [272]) Assumptions and notation of Problem \2.2 in force. 
Consider the stochastic differential equation 

dYt = uUYt)dt + at{Yt)dWt, lo ~ Po, 



(3) 



dt 



d'] r (c(0 - Vem)) exp[e'MO - c{x)M, 

/ J — oo 



^In order to contain space and notation the underlying geometric setup is not fully developed here. We just 
say that the problem originated from the use of differential geometry and statistics for the nonlinear filtering 
problem. The reader interested in geometric aspects and other details is referred to Brigo (1996), Brigo, Hanzon 
and Le Gland (1999), or to the tutorial in Brigo (1999). 
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where V0ip{9t) = {dip / 09^ (Ot) , . . . , dip / dO^{6t)Y , with the symbol "~" to he read as "distributed 
as". 

Ifu" e U{pq,(j), then the SDE ^ solves Problem{KM in that 

PYAx) = exp[e'tc{x)-ip{et)], t>o. 

The proof of the theorem is rather straightforward. It is sufficient to write the Fokker-Planck 
equation for the SDE ([3]) and, after lengthy computations, verify that indeed 

- exp[^^c(x) - m)] = exp[e[c{x) - m)]) + i^ti^) exp[e'Ax) - m)]) 

by substituting the expression for u given in the theorem. A different proof can be found in 
Chapter 7 of Brigo (1996) or in Brigo (2000), where in deriving the expression for u it was 
tacitly assumed, as is done here, that 

lim u'p{x)pt{x) = for all t > 0. 



x— >— oo 



In the next section, we shall consider an interesting application of this theorem to the option 
pricing problem. Indeed, we shall use such result more as a "guiding tool" rather than applying 
it immediately as it stands. In particular, assumptions (Al) and (A2) will be checked directly 
and not via the sufficient conditions usually considered in the literature. 

3 Alternatives to the Black and Scholes model 

Let us consider the Black and Scholes (1973) stock price model, 

dSt = ^iStdt + aStdWt, So = so, tG[0,T], (4) 

where sq is a positive deterministic initial condition, and fi, a and T are positive real constants. 
The probabihty density pst of St, at any time t > 0, is given by 



PSt[X) 



exp|cin- + p(t)ln2^-^(C,p(t))|, a; > 0, (5) 



^2 2' ' 2aH' 

With the notation for exponential families introduced in the previous section, one writes 

ci(a;) = In -, c^^x) = In^ -, Ot = {C, p(t)}' . 

So So 
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One might wish to model the stock price process by considering a different volatility function 



cr, instead of a St in (|l]o, while preserving major properties of the original process The 
purpose of this section is then the construction of alternative stock price dynamics that differs 
from (jlj), yet sharing similar features from a probabilistic point of view. 

Let us approach this problem by applying Theorem 12.31 to find a SDE with a given diffusion 
coefficient (Jt{-) and whose marginal density is equal to the marginal density of S in all time 
instants of the time interval T = [e, T], where < e < T and e can be chosen arbitrarily close to 
0. We remark that the density of S is not concentrated in the whole real line, but in (0, +oo). 
However, the procedure resulting in Theorem 12.31 can be adapted in a straightforward way to 
the latter case once —oo is replaced by 0. 

It easily follows from ([3]) that the equation sought for is 



dYt = uUYt,So,0)dt + at{Yt)dWt, Y, = S,, e < t < T, 

X C + 1 



(6) 



C + 2p(t-a)ln- 

y 



X 

+ 



2{t-a) 



In 



y 2p{t-a) 



dx X 

The definition of Y is then extended to the whole interval [0, T] by setting 

dYt = jjYt dt + aYtdWt, < t < e, Yq = sq . 

In other terms, Y is assumed to follow the same dynamics as the original process S in [0,e), 
and to start from the same initial condition. 

The reason for splitting the interval [0,T] into two subintervals is simply to avoid problems 
concerning the definition of the drift at time t = 0. The variable a is here introduced to allow 
for subsequent generalizations. 

In general, we cannot prove that ([6]) satisfies assumptions (Al) and (A2) for a prescribed 
cr, i.e. that G U{pQ,a) so that lA{pQ,a) is nonempty. Therefore, at this stage, ^ is simply 
providing us with a candidate solution of our problem for a generic a. 

In the following table we report some possible choices for a and the corresponding candidate 
. In particular, we can notice that for (Tj(x) = ax, So,0) equals, as it must be, the 

original drift px of the Black and Scholes model (jlj). 

Let us now assume that the particular a we are working with is such that u" G U{pQ, a). 

We have seen so far that, by applying the results of the previous section, it is rather straight- 
forward to produce processes with the same marginal distribution as that of (jlj). However, a 
further fundamental property of process (jlj) is that its log-returns between any two time instants 
are normally distributed, independently of the prices at the considered instants, i.e., 

^ M Up - \(j^)5, (jH) , 5>0, tG[0,T-5]. (7) 



^We use the term "volatility" to denote the whole diffusion coefficient crt(-) rather than the standard deviation 
rate of the instantaneous return as usually done in practice. 
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Table 1 



Examples of volatility functions and corresponding candidate drift 





Drift 7;°'('r Sn 


V 


J C + 2p(t) In ^ + I In ^ ,^+\ 




4(^-1) -^f^ln^ + f(/i-f) + iln^ 


vx 


^ i(^^-^^) + f(S + l) +iln|;(l-S) 


ax 


\ix 



Alternative models such as ([6]) do not share this property in general. In fact, identity of the 
marginal laws alone does not suffice to ensure ([7]) , for which equality of second order laws or of 
transition densities would be sufficient instead. How can we obtain alternative models whose 
properties concerning log- returns are as close as possible to property ([7])? 

To tackle this issue, we have to find a compromise between our alternative model ([6]) and 
model (jl]). To this end, we consider a weaker version of ([7]) in that we restrict the set of dates 
for which the property holds true. Precisely, we modify the definition of Y so that, given the 
time instants := {0, A, 2A, . . . , A^A}, A = TjN, A > e, property (JTl) is satisfied by Y in 
r^, i.e. 

ln^~Ar((p-ia2)(^-j)A, a2(2-j)A), i> 3. ^ = l,...,iV, j = 0, . . . , AT - 1. (8) 

Limiting such key property to a finite set of times is not so dramatic. Indeed, only discrete time 
samples are observed in practice, so that once the time instants are fixed, our process Y can 
not be distinguished from Black and Scholes process'. The issue of discrete versus continuous 
time will be further developed in Section 5. 

The new definition of Y is still based on Theorem 12.31 However, we use this theorem 
"locally" in each time interval [(i — 1)A, iA). This means that in such interval we define 
iteratively the drift u'^ as in the theorem but 

• we translate back the time-dependence of a time amount (i — 1) A (thus locally restoring 
the dynamics of the original result) and 

• we replace the distribution for the initial condition with the distribution of the final 
value of Y relative to the previous interval. 

We obtain: 

dYt = u'l{Yt,Y,^t),ait))dt + atiYt)dWt, t e [lA + e, (t + 1)A), (9) 

dYt = fiYtdt + aYtdWt, for te[iA,tA + e), a{t) = lA for te[iA, (z + 1)A) , 
where u^{x,y,a) was defined in ([6]). 
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It is clear by construction that the transition densities of S and Y satisfy PY^i+i^^\YiAi^i v) = 
PS(i+i)A\SiAi^'j y)- Then, starting from the equahty of the marginal laws of S and Y in the first 
interval that holds by construction, we inductively obtain the equality of the marginal laws 
also in each other interval. As a consequence, the second order densities are also equal among 
consecutive instants {i — 1)A, -iA, i.e.. 

It follows that 

ln%ll^~Ar((/.-|a>, a^A), z = 0, . . . , iV - 1. 

At this point we remark that the process y in ([9]) is not a Markov process in [0, T] 
it is Markov in all time instants of T^. Formally, 

Y^m-2)A,-,Yo - PY^a\Y( 

This property follows from the fact that in [(m — l)A,mA) the dynamics of the SDE defining 
Y does not depend on Y(m-2)A, • • • , ^o, and that when such equation is considered for t G 
[(m — l)A,mA), in its drift u" the local initial condition for the entry Y is set to y(m-i)A- 
From now on, we refer to markovianity in as to A-Markovianity. 

We finally notice that, through the A-Markovianity, property (fTOj) extends to any pair of 
instants in T^, so as to yield (IHl). Moreover, the inductive application of the A-Markovianity 
and the identity of transition densities in the grid leads to the identity of the finite dimentional 
distributions of S and Y in the grid. 

4 Option pricing in continuous-time 

Let us now consider the process {Bt : t > 0} whose value evolves according to 

dBt = Btrdt, (11) 

with Bq = 1 and where r is a positive real number, so that Bt = exp(rt). The process B is 
assumed to describe the evolution of a money market account in a given financial market. The 
process Yt in Qj is instead assumed to model the evolution of some traded financial (risky) 
asset, typically a stock. 

The financial market thus defined might admit arbitrage opportunities. As is well known, 
a sufficient condition which ensures arbitrage-free dynamics is the existence of an equivalent 
martingale measure with respect to the initially chosen numeraire. In this paper, we use the 
process i? as a numeraire, so that an equivalent martingale measure is a probability measure 
that is equivalent to the initial one, P, and under which the process {Yt/Bt : t > 0} is a 
martingale. A necessary condition for the existence of an equivalent martingale measure is the 
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semimartingale property for the process Y. The process Y is indeed a semimartingale under P 
for sufficiently well behaved volatility functions crj(-). 

We denote by S the set of all volatility functions at{-) such that u'^ G U{po, a) and for which 
there exists a unique equivalent martingale measure. 

The set S is obviously non-empty, since it contains at least the Black and Scholes volatility 
function at{x) = ax. Moreover, as we will prove in the sequel, all volatilities functions of type 
i>I, u > 0, belong to S, with / denoting the identity map. An interesting example of volatility 
functions which do not belong to S is instead provided in the appendix. 

We now assume that we have chosen a G iS and the corresponding equivalent martingale 
measure Q'^. Since {Yt/Bt : t > 0} is a martingale under such a measure, it easily follows that 
under the process Y satisfies the SDE 

dYt = rYt dt + aYt dWt, t G [iA, iA + e), 

dYt = rYt dt + at{Yt) dWt, t G [zA + e, (z + 1)A), 

where is a standard Brownian motion under Q'^. 

Furthermore, under the assumption that i) there are no-transaction costs, ii) the borrowing 
and lending rates are both equal to r, iii) short selling is allowed with no restriction or penalty, 
and iv) the stock is infinitely divisible and pays no dividends, the unique no-arbitrage price for 
a given contingent claim H G L'^{Q") is (see Harrison and Pliska (1981, 1983)) 

Vt = ^E<^^ {H\Tt}^ (12) 

where {J-'t : t > 0} denotes the filtration associated to the process Y . 

In the special case of a European call option, the following are interesting problems to solve. 



Problems 4.1 Let us assume that the given claim is a European call option, written on the 
stock, with maturity T and strike K . Find: 



■mi B^'E^^ [{Yt-K)^\:F,}., 
sup B^^E^^ I {Yt - K)+ Joj 



(13) 
(14) 

e>0,aeS 

Solving these problems is equivalent to finding the lowest and highest theoretical price of the 
option for which the underlying stock price has lognormal marginal distribution and normal 
log- returns on the grid T^, with standard deviations proportional to a. 

If we denote by K the value of the infimum in (fT3|) . and by V* the value of the supremum 
in (fT4l) . the following inequalities obviously hold: 



(sq - Ke 



< K < VBs{(r) <V*< So, 



(15) 
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where Vbs((7') denotes the option Black and Scholes price at time as determined by a, the 
volatihty parameter in (jl]). Indeed, since a I G S, the central inequalities hold by definition of 

and V*, whereas the first and the last ones feature respectively the well known no-arbitrage 
lower and upper bounds for option prices. 

In the next subsection we will show that the first and last inequalities in (ITSjl are actually 
equalities. To prove this statement, it will be sufficient to restrict our analysis to the class of 
volatilities {crf(x) = ux, v > 0}. This result is at first sight surprising. Indeed, one would 
naively expect that the difference between prices implied by models which are equivalent in 
the A-grid is bounded by a quantity that is somehow related to A, typically (9(A'*') for some 
positive real A. In fact, by halving the size of A, we double the discrete-time instants where 
the models in our family are equivalent. Accordingly, we would expect the prices implied by 
these now "closer" models to range in a narrower interval. However, as we shall soon see, this 
is not the case. 

4.1 A fundamental case 

We begin by stating the following. 

Lemma 4.2 In the fundamental case where (Jt{x) = vx, > 0, the process given by 
dY^" = uf{Yt, Y^^t),a{t))dt + iyYt dWt, t e [lA + e, (t + 1)A), 



(16) 



uf{y,ya,a) = y 



U 2 -2n , ^/'^^ 



■(l-^)ln:^^, 



2(t-a) cr2 y^' 

dYI" = jj^Y^^dt + aY^^dWi for te[tA,iA + e), a{t) = iA for te[iA, {i + l)A) 

solves Problem\EMwhen pt is given by Moreover, the volatility function <Jt{x) = vx belongs 
to S, for any u > 

Proof. Since Y^ has the same marginal distribution as S under P, it follows that Y^ > 0. Then, 
the process Zt = InY^" is well defined and, by Ito's formula, the SDE for Zt is piecewise linear 
in a narrow sense (in that its diffusion coefficient is purely deterministic), and hence admits a 
unique strong solution which, for each t G [j'A, (j + 1)A), is explicitly given by 

Zt = Z,A + (/x-ia2)(t-jA) (17) 
a{Wt-W,A) te[jA,jA + e), 



t-jA 



[3/2 



te[jA + e, (j + l)A), 



,,2 



where P = 1 — jj- 

As a consequence, the assumptions of Problem l2.2l are satisfied so that u'^^ solves Problem l2.2l 
when pt is given by ([5]) . Moreover, the Girsanov change of measure from P to Q"^ is well defined 
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since one can show that the Novikov condition is satisfied through apphcation of the "tower 
property" of conditional expectations. Hence, the measure Q'^^ exists unique and ul belongs 
to S for each u > 0. For the uniqueness of the measure Q"^ , we refer for instance to Duffie 
(1996). 

Theorem 4.3 In the fundamental case where crt(x) = ux, u > 0, the unique no-arbitrage 
option price at time t is the Black and Scholes price 



where 



U\t, v) = Y.^^d,) - Ke-'^^-''><!>{d2), 



^ UXr/K) + (r + z7^(t)V2)(T - t) 



^it)VT^t 



d2 = di{t) - v\t)^T -t, 

u^t) = { \ ^ (19) 

^ lEZEE^pE^ tGKt) + e,a(t) + A) 

Proof. From the previous lemma, we infer the existence of a unique no-arbitrage option price 
that can be calculated through f|T2l) . Then f|T8l) is obtained by noticing that under the equivalent 
martingale measure 

^ M {{r - \u\tf){T - t\ z7^(t)2(T-t)), tG[0,T], 

with z/"^ given by (|T9l) . This implies that the option price at time t corresponding to is the 
Black and Scholes price with volatility coefficient z?'^(t), i.e., that ( JTSll holds. 



The key point of our result is that, for any given volatility coefficient z/, we are free to adjust the 
drift of the SDK defining the dynamics of the stock-price process under the objective measure, in 
such a way that the resulting Y^ has the same distributional properties of the Black and Scholes 
process on discrete-time dates. As opposed to this, the risk-neutral valuation for pricing options 
imposes the drift rY" to the SDK followed by Y^ under the equivalent martingale measure. This 
causes the option price implied by the alternative model to coincide, at first order in e, with 
the Black and Scholes price with volatility parameter v. In fact, imposing the drift rY'^ to Y'^ 
leads to the same risk neutral process as that of Black and Scholes' (obtained by imposing the 
drift rS to S*), with the only difference that a is replaced by v. 

Remark 4.4 (Historical versus Implied volatility). A first interesting property that can 
he deduced from this theorem applies in case one believes that option prices trade independently 
of the underlying stock price. We have in fact been able to construct a stock price process, 
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the process fT^) . whose marginal distribution and transition density depend on the volatility 
coefficient a, whereas the corresponding option price, in the limit for e 0, only depends on 
the volatility coefficient v. As a consequence, we can provide a consistent theoretical frame- 
work which justifies the differences between historical and implied volatility that are commonly 
observed in real markets. 

Straightforward application of the previous theorem leads to the main result of this section 
which is summarized in the following. 

Corollary 4.5 The solutions of problems [TB^] and are 

K = (so - Ke-^'^y 

(20) 

V* = So. 

Moreover, for any other candidate price V G {V^,V*) there exist a volatility v and an e > 
such that 

U'{0,iy) = V. 

Proof. To prove fl20|) . we simply have to take the limit of expression f fTSj) (with t = 0) for e 
going to zero and u either going to zero or going to infinity, since 



\im Vbs{v) = {so - Ke 



lim Vbs{v) = So, 

and 

lim z/^(0) = u. 

Finally, we remember that, ceteris paribus, VBsip) is a strictly increasing function of f , so 
that VBsip) = y has a unique solution for V G (K, V*), hence t/^(0, u) = V has a solution in 
(0,+oo) X {V,,V*). 

Remark 4.6 (Taking e ^ ). It is possible to consider the limit for e —* in the above 
expressions so as to present our result in a simpler and more elegant way. However, the 
treatment with e does not involve limit considerations and permits to contain analytical effort, 
so that we decided to keep e > 0. This can be also useful in numerical implementations. We 
just observe that, for e 0, since (3 < 1, becomes 



Zt = Z,A + (/i-ia^)(t-jA) 

+(t - jAf/'u f {u - jA)-P'^dWu t G [jA, (j + 1)A). 

This process is well defined since the integral in the right-hand side exists finite a.s. even though 
its integrand diverges when u j A^. 
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The above equation can be better compared to the Black and Scholes process when written 
in differential form: 

dZt = (/i - \a^) dt + f^{t- jA)^/2-V /* {u - 3/\)-^'^dWu dt + udWt te [jA, (j + 1)A). 

Z JjA 

By observing this last equation we can isolate three terms in the right-hand side. The first term 
is the same drift as in the log-returns of the Black and Scholes process. The third term is the 
same as in the log-returns of the Black and Scholes process, but the volatility parameter a is 
replaced by our v. Finally, the central term is the term which is needed to have returns equal 
to the returns in the Black and Scholes process even after changing the volatility from o to v. 
Note that this term goes to zero for a = v. It is this term that makes our process non-Markov 
outside the trading time grid. 

The interpretation of the previous theorem and corollary is as follows. If we are given a 
discretely observed stock price, the particular way we use to "complete" the model with any 
of our continuous-time processes Y has a heavy impact on the associated option price. The 
influence is in fact so relevant that such a price can be arbitrarily close to either no-arbitrage 
bounds for option prices. 

From an intuitive point of view, the reason why option prices can be so different is because 
the time step in the grid is never infinitesimal. In other words, continuous-time option 
prices would simply reveal the differences existing at an infinitesimal level among all the stock 
price processes . 

5 Option pricing and hedging in the real world 

Let us now consider a trader who needs to price an option on a given stock. His usual practice is 
to resort to continuous-time mathematics to fully exploit the richness of its theoretical results, 
and to model stock returns with a normal distribution. 

Let us denote by 5t the length of the smallest time interval when an actual transaction can 
occur. Such 5t is the best realistic approximation of the infinitesimal time distance "cit" . 

The results of the previous section imply that the geometric Brownian motion (jlj) is just 
one of the infinitely many processes Y that possess the properties required by the practitioner 
along intervals of equal length 5t. However, the basic equivalence in the description of the 
stock price dynamics can not be extended to the corresponding option prices. Indeed, any real 
number in the interval (K, V*) can be viewed as the unique no-arbitrage option price for some 
process Y . 

Which process Y should then be chosen by the practitioner? 

A possible answer to this question can be provided, for example, through the estimation 
of the option replication error in discrete time. To this end, for any z/ > 0, let us denote by 
{i^iV'^) the self-financing strategy that replicates in continuos time the option payoff for the 
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process V^, where and rj^ are respectively interpreted as the number of stock shares and 
money units held at time t. 

Then, fixing a set of dates tq = < ri < ■ • ■ < r„ = T such that = iSt% and, denoting 
the observed stock price at time tj by Sr^, j = 0,. . . ,n, the replication error when hedging 
according to the strategy {C,'^,i]'^), starting from the endowment t/'^(0, z/), is 

n— 1 n—1 
j=0 j=0 

At this stage, we can solve, for example, an optimization problem where e(z/) is minimized, 
according to some criterion, over all u > 0. Such procedure, however, is justified only to measure 
the performance of our continuous-time prices and strategies on the fixed set of discrete-time 
instants. More generally, the issue of deriving a fair 6t-time option price and hedging strategy 
should be tackled by resorting to the existing literature on incomplete markets. 

Many are the criterions one can choose from for pricing and hedging in incomplete mar- 
kets. We mention for instance those of FoUmer and Sondermann (1986), FoUmer and Schweizer 
(1991), Schweizer (1988, 1991, 1993, 1994, 1995, 1996), Schal (1994), Bouleau and Lamberton 
(1989), Barron and Jensen (1990), El Karoui, Jeanblanc-Picque and Viswanathan (1991), Davis 
(1994), Frasson and Runggaldier (1997), El Karoui and Quenez (1995), FritteUi (1996), Mer- 
curio (1996), Mercurio and Vorst (1997), Bellini and FritteUi (1997), Frey (1998) and FoUmer 
and Leukert (1998). 

However, the purpose of this section is not to favor any particular approach. We want, 
instead, to stress the following innovative feature in the theoretical problem of option pricing. 
Instead of fixing a stock price process and then deriving a fair option price and an "optimal" 
hedging strategy, we can in fact consider a family of processes, that are somehow equivalent 
in the description of the stock price evolution, among which we can select a convenient one by 
means of our favorite incomplete markets criterion. 

A comparison between the performances of the strategy (^'^, t]") and the hedging strategy 
associated to any incomplete-market criterion is beyond the scope of the paper and is left to 
future research. 

6 Conclusions 

In the present paper we consider an option-pricing application of a result which is based on the 
construction of nonlinear SDE's with densities evolving in a given finite-dimensional exponential 
family. Precisely, we derive a family of stock price models that behave almost equivalently to 
that of Black and Scholes. All such models share the same distributions for the stock price 



^ Although fixed a priori, the set T'^ introduced earfier is arbitrarily chosen, so that we can set A — St/k, 
with k any positive integer implying that {ti, r„} C T*^. 



D. Brigo, F. Mercurio. Banca IMI - PDG internal report 



16 



process and its log- returns along any previously fixed 'trading time-grid'. Therefore, all these 
models can be viewed as equivalent in the description of the stock price evolution. 

However, the continuous-time dynamics chosen to 'complete' the model in between the 
instants of the trading time-grid, reflects heavily on the option price. The option price in fact 
can assume any value between the option intrinsic value and the underlying stock price. 

As a conclusion, our result points out that no dynamics is the right one a priori, and that an 
incomplete-market criterion is needed to choose among all the different models. Practitioners 
with different criteria can still pick up a model from our family, so as to match their expectations 
or to minimize their exposures. 

Even though our results are based on the assumption of a lognormal distribution for the stock 
price and a normal distribution for its log-returns, the generalization to many other distributions 
is possible. Notice, indeed, that the results of Section 2 can be applied to any curve of densities 
in a given exponential family. However, the generalization is not straightforward and heavily 
relies on the particular distributions which are considered. 



Appendix 



In this appendix we consider the case crt{') = 7^ as an interesting example of a class of 
volatilities which do not belong to S. Suppose, for the sake of simplicity, that we require only 
the marginal law of Y to coincide with the marginal law of S, ignoring the returns distributions. 
Now assume by contradiction that v & S. Then by applying ([6]) we obtain a Markov process Y 



dY. ^ \- 



C + 2p{t) log ^ 



dt+ — 
2t 



1 ^* 
log — 

So 



C + 1 

2p{t) 



dt + u dWt , e<t<T , 



(21) 



with Yt = St for each t G [0, e]. 

Let us focus on t > e. Since we now know that Yt and St have the same distribution under 
the objective probability measure, Yt is lognormally distributed, and in particular Yt > for all 
t. Since z/ G 5, there exists an equivalent martingale measure, so that, under such a measure, 

dYt = rYt dt + u dWt , 

where is a standard Brownian motion under the martingale measure. Now notice that this 
last SDE is a linear equation, so that its solution has a density whose support is the whole 
real line. In other terms, such solution can be negative with positive probability at any fixed 
time instant, as opposed to what we have seen for Y under the objective measure. However, 
there cannot be an equivalent measure that transforms a process whose support is the positive 
halfline into another one whose support is the whole real line. Therefore, we have contradicted 
our assumption that z/ G 5. 

Our example is similar in spirit to that of Delbaen and Schachermayer (1995) who consider a 
Bessel process (taking positive values at all times) which cannot be transformed into a Brownian 
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motion. Some conditions for the existence of an equivalent martingale measure are given for 
example in Rydberg (1997). 

References 

[1] BARRON, E. N. and JANSEN, R. (1990). A Stochastic Control Approach to the Pricing 
of Options. Mathematics of Operations Research 15, 49-79. 

[2] BLACK, F. and SCHOLES, M. (1973). The Pricing of Options and Corporate Liabilities. 
Journal of Political Economy 81, 637-654. 

[3] BELLINI, F. and FRITTELLI, M. (1997). Certainty Equivalent and No Arbitrage: a 
Reconciliation via Duahty Theory. Working Paper. Brescia University. Italy. 

[4] BOULEAU, N. and LAMBERTON, D. (1989). Residual Risks and Hedging Strategies in 
Markovian Markets. Stochastic Processes and their Applications 33, 131-150. 

[5] BRIGO , D. (1996). Filtering by Projection on the Manifold of Exponential Densities. PhD 
Thesis, Free University of Amsterdam. 

[6] BRIGO, D. (1997). On nonlinear SDEs whose densities evolve in a finite-dimensional 
family, in: I. Csiszar and Gy. Michaletzky (Editors), Stochastic Differential and Difference 
Equations, Birkhauser, Boston, 11-19. 

[7] BRIGO, D. (1999) Diffusion Processes, Manifolds of Exponential Densities, and Nonlin- 
ear Filtering, in: O.E. Barndorff-Nielsen and E. B. Vedel Jensen (Editors), Geometry in 
Present Day Science, World Scientific. 

[8] BRIGO, D. (2000). On SDEs with marginal laws evolving in finite-dimensional exponential 
families. Statistics and Probability Letters 49 , pp. 127-134 

[9] BRIGO, D., HANZON, F., and LE GLAND, F. (1999). Approximate nonlinear filtering 
by projection on exponential manifolds of densities. Bernoulli, Vol. 5, N. 3, pp. 495-534. 

[10] BRIGO, D., and MERCURIO, F. (2000). Option Pricing Impact of Alternative 
Continuous-Time Dynamics for Discretely-Observed Stock Prices. Finance and Stochas- 
tics. Vol. 4, N. 2. 

[11] DAVIS, M. H. A. (1994). A general Option Pricing Formula. Preprint. Imperial College, 
London. 

[12] DELBAEN, F., and SCHACHERMAYER, W. (1995). Arbitrage possibilities in Bessel 
processes and their relations to local martingales. Probability Theory and Related Fields 
102, 357-366. 



D. Brigo, F. Mercurio. Banca IMI - PDG internal report 



18 



[13] DUFFIE, D. (1996). Dynamic Asset Pricing Theory, Second Edition. Academic press, New- 
York. 

[14] EL KAROUI, N., JEANBLANC-PICQUE, M. and R. VISWANATHAN (1991). On tlie 
Robustness of Black-Scholes equation. Laboratoire de Mathematiques et Modelisations, 
Ecole Nationale des Fonts and Chaussees-Ecole Normale Superieure de Cachan, Prepubli- 
cation 91/1. 

[15] EL KAROUI, N., and QUENEZ, M. C. (1995). Dynamic Programming and Pricing of 
Contingent Claims in an Incomplete Market. SIAM Journal on Control and Optimization 
33, 29-66. 

[16] FOELLMER, H. and SONDERMANN, D. (1986). Hedging of Non-Redundant Contingent 
Claims. In W. Hildenbrand and A. Mas-Colell (eds.). Contribution to Mathematical Eco- 
nomics, 205-223. 

[17] FOELLMER, H. and SCHWEIZER, M. (1991). Hedging of Contingent Claims Under 
Incomplete Information. In Applied Stochastic Analysis (M.H.A. Davis and R.J.Elliott 
eds.). Stochastic Monographs vol.5, Gordon & Breach, London/New York, 389-414. 

[18] FOELLMER, H. and LEUKERT, P. (1998). Quantile Hedging. Preprint. Humboldt- 
Universitat Berlin. 

[19] FRASSON, S. and RUNGGALDIER, W.J. (1997). A Stochastic Control Model for Hedging 
in Incomplete Markets. Statistics and Control of Stochastic Processes, edited by Yu. M. 
Kahanov, B.L.Rozovskii and A.N. Shiryaev. World Scientific. Singapore. 

[20] FREY, R (1998). Superrephcation in Stochastic Volatility Models and Optimal Stopping. 
Preprint. Department of Mathematics. ETH Ziirich. Switzerland. 

[21] FRIEDMAN, A. (1975). Stochastic Differential Equations and Applications, I. Academic 
Press, New York. 

[22] FRITTELLI, M. (1996). The Minimal Entropy Martingale Measure and the Valuation 
Problem in Incomplete Markets. Journees Internationales de Finance, AFFI, Geneva, 
June 1996 

[23] HARRISON, J.M. , and PLISKA, S. R. (1981). Martingales and Stochastic Integrals in the 
Theory of Continuous Trading. Stochastic Processes and Their Applications, 11, 215-260. 

[24] HARRISON, J.M. , and PLISKA, S. R. (1983). A Stochastic Calculus Model of Continuous 
Trading: Complete Markets. Stochastic Processes and Their Applications, 15, 313-316. 

[25] MERCURIO, F. (1996). Claim pricing and hedging under market imperfections. PhD 
Thesis, Tinbergen Institute, Rotterdam. 



D. Brigo, F. Mercurio. Banca IMI - PDG internal report 19 

[26] MERCURIO, F. and VORST A.C.F. (1997), Option Pricing and Hedging in Discrete Time 
with Transaction Costs. Mathematics of Derivative Securities, edited by M.A.H. Dempster 
and S.R. Pliska, 190-215. 

[27] MERTON, R. (1973). Theory of Rational Option Pricing. Bell Journal of Economics and 
Management Science 4, 141-183. 

[28] ROGERS, L.C.G. and WILLIAMS, D. (1987). Diffusions, Markov Processes, and Martin- 
gales, VoL II. Wiley & Sons. 

[29] ROGERS, L.C.G. and SATCHELL, S.E. (1996). Does the Behaviour of the Asset Tell us 
Anything about the Option Price Formula? A Cautionary Tale. Preprint. University of 
Bath. 

[30] RYDBERG, T. H. (1997). A note on the existence of unique equivalent martingale measures 
in a Markovian setting. Finance and Stochastics 4, 251-257. 

[31] SCHAL, M. (1994). On Quadratic Cost Criteria for Option Hedging. Mathematics of Op- 
eration Research 19, 121-131. 

[32] SCHWEIZER, M. (1988). Hedging of Options in a General Semimartingale Model. Diss. 
ETHZ No 8615. 

[33] SCHWEIZER, M. (1991). Option Hedging for Semimartingales. Stochastic Processes and 
their Applications 37, 339-363. 

[34] SCHWEIZER, M. (1993). Approximating Random Variables by Stochastic Integrals, and 
Applications in Financial Mathematics. Habilitationsschrift. University of Gottingen. 

[35] SCHWEIZER, M. (1994). Risk-Minimizing Hedging Strategies Under Restricted Informa- 
tion. Mathematical Finance 4, 327-342. 

[36] SCHWEIZER, M. (1995). Variance Optimal Hedging in Discrete Time. Mathematics of 
Operation Research 20, 1-32. 

[37] SCHWEIZER, M. (1996). Approximation Pricing and the Variance Optimal Martingale 
Measure. Annals of Applied Probability 24, 206-236. 

[38] STROOCK, D.W. and VARADHAN, S.R.S (1979). Multidimensional diffusion processes. 
Springer, Berlin. 



